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1^ (57) Abstract: Multi-slice magnetic resonance imaging of a region of interest of an imaging subject (16) is performed using a radio 
frequency coil (40) arranged to generate a Bi magnetic field in the region of interest. One or more processors (44, 82, 88, 110) 
determine a Bl field value for each slice that is representative of the Bl field over a selected area of the slice, accounting for subject 
effects on the BI field, and determine an adjusted per- slice radio frequency excitation for each slice that adjusts the Bi field value for 
the slice to a selected value. A magnetic resonance imaging system (10, 44, 46, 50, 52) acquires magnetic resonance imaging data 
for each slice using the adjusted per-slice radio frequency excitation for that slice. A reconstruction processor (58) reconstructs the 
acquired magnetic resonance imaging data into a reconstructed image representation. 
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1 

Bl FIELD CONTROL IN MAGNETIC RESONANCE IMAGING 

DESCRIPTION 

The following relates to the magnetic resonance arts. It finds particular 
application in magnetic resonance imaging, and will be described with particular reference 
thereto. However, it also finds application in other magnetic resonance applications such as 
magnetic resonance spectroscopy. 

5 As the flip angle of spins relates to the Bl magnetic field, it is desirable that 

the Bl magnetic field produced by radio frequency transmissions during magnetic 
resonance imaging be spatially uniform, both within an image slice, and also across 
multiple image slices. At relatively low resonance firequencies such as 64 MHz 
corresponding to a main (Bo) magnetic field of about 1.5T, modifications of the Bi field by 

0 conductive and dielectric characteristics of the imaging subject are typically negligible. 
However, for firequencies of 128 MHz or higher (Bo > -3.0T for ^H), conductive and 
dielectric characteristics of the imaging subject substantially affect the Bi field uniformity 
and cause the field to be more non-uniform than it would be without subject effects. 

To partly mitigate the effects of Bi field non-uniformity, an adiabatic radio 

5 firequency pulse packet can be employed. The adiabatic radio firequency pulse packet 
reduces the effect of spatial variations in radio frequency field strength to produce a more 
uniform flip angle. However, adiabatic radio frequency pulse packets have certain 
disadvantages, including being difficult to refocus and employing a hi^er radio frequency 
power levels resulting in higher specific absorption rate (SAR) values. 

0 An approach for partly mitigating the effects of Bi field non-uniformity 

across multiple slices is disclosed in Clare et al., "Compensating for Bi inhomogeneity 
using active transmit power modulation," Magnetic Resonance Imaging vol. 19, pp. 
1349-52 (2001). To improve uniformity of flip angles across multiple axial slices, Clare et 
al. measured the axial Bi field dependence along the central axis of an empty birdcage 

5 radio frequency coil using a magnetic flux probe to produce a lookup table of radio 
frequency excitation intensity correction values. Using the lookup table, the radio 
frequency power input was adjusted on a per-slice basis for each acquired slice to 
compensate for the axial dependence of the Bi field. Measurements at 128 MHz (3.0T) for 
homogeneous gel-filled cylindrical phantoms and human heads using this input power 
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compensation produced a substantial axial uniformity improvement in flip angle along the 
central axis of the coil. However, substantial in-slice non-uniformity remained, and it was 
recognized that the axial correction along the central axis is unable to compensate for such 
in-plane Bi field variation. 
5 At higher frequencies, ttie Bi fields in the human body are known to differ 

substantially from those measured in an empty or lightly loaded RF coil. Moreover, 
because the Bi field variation across slices is expected to depend upon the in-slice 
coordinate and effects of the subject, a calibration generated along the central axis of an 
empty coil may be largely ineffective for compensating variation of the Bi field at radial 
10 positions away from the central axis along which the axial flux probe calibration was 
acquired. 

The present invention contemplates an improved apparatus and method that 
overcomes the aforementioned linoitations and others. 

15 

According to one aspect, a magnetic resonance imaging method is provided 
for performing multi-slice magnetic resonance imaging of a region of interest of an 
imaging subject using a radio frequency coil arranged to generate a Bi magnetic field in 
the region of interest. A per-slice Bi field value is determined for each slice that is 

20 representative of the Bi field over a selected area of the slice. An adjusted per-slice radio 
frequency excitation is determined for each sUce that adjusts the Bi field value for the slice 
to a selected value. Magnetic resonance imaging data are acquired for each slice using the 
adjusted per-slice radio frequency excitation for that slice. The acquired magnetic 
resonance imaging data are reconstructed into a reconstructed image representation. 

25 According to another aspect, a magnetic resonance imaging apparatus is 

disclosed for performing multi-slice magnetic resonance imaging of a region of interest of 
an imaging subject. A radio frequency coil is arranged to generate a Bi magnetic field in 
the region of interest. A means is provided for determining a per-slice Bi field value for 
each slice that is representative of the Bi field over a selected area of the slice. A means is 

3 0 provided for determining an adjusted per-slice radio frequency excitation for each slice that 
adjusts the Bi field value for the slice to a selected value. A means is provided for 
acquiring magnetic resonance imaging data for each slice using the adjusted per-sUce radio 
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frequency excitation for that slice. A means is provided for reconstructing the acquired 
magnetic resonance imaging data into a reconstructed image representation. 

One advantage resides in providing a radio frequency power input 
correction on a per-slice basis that accounts for in-slice variation of the Bi field. 
5 Another advantage resides in providing Bi field correction that accounts for 

spatially varjdng conductance and dielectric properties of the imaging subject. 

Another advantage resides in providing Bi field correction in combination 
with determination of peak and average specific absorption rate (SAR) values to optimize 
the Bi field for both uniformity and intensity while ensuring the SAR remains within 
1 0 regulatory limits. 

Yet another advantage resides in providing a per-slice Bi field correction 
based on a simulation of the Bi field over an area of interest of the slice, where the 
simulation is performed using a nonhomogeneous model of at least a region of interest of 
the imaging subject. 

15 Still yet another advantage resides in calibrating the radio frequency power 

input based on in situ measurement of the Bi field over an area of interest of the slice with 
the imaging subject or with a compartmentalized non-homogeneous phantom load. 

Numerous additional advantages and benefits will become apparent to those 
of ordinary skill in the art upon reading the following detailed description of the preferred 

2 0 embodiments. 



The invention may take form in various components and arrangements of 
components, and in various process operations and arrangements of process operations. 
25 The drawings are only for the purpose of illustrating preferred embodiments and are not to 
be construed as limiting the invention. 

FIGURE 1 shows a magnetic resonance imaging system including a 
per-slice calibration of radio frequency excitation. 

FIGURES 2A and 2B show a top view and a side section, respectively of 
30 the TEM head coil of FIGURE 1. In the top view of FIGURE 2 A, the rungs are shown in 
phantom. The side section of FIGURE 2B is taken along a plane passing through two 
maximally spaced-apart rungs. 
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FIGURES 3A, 3B, and 3C show axial, sagittal, and coronal slices, 
respectively, of a TEM head coil model with a human head model disposed therein, for 
performing electromagnetic simulation of the Bi field produced by the TEM head coil. 

FIGURE 4 plots the average |Bi| field across a slice versus axial position z 
5 obtained using electromagnetic simulation of the model shown in FIGURES 3A, 3B, and 
3C, and using electromagnetic simulation using the same head size but using averaged 
tissue properties. 

FIGURE 5 plots normalized \Bi\ field standard deviation in slice versus 
axial position z for the realistic head model of FIGURES 3 A, 3B, and 3C, and for a head 
1 0 model with averaged tissue properties. 



With reference to FIGURE 1, a magnetic resonance imaging scanner 10 
includes a housing 12 defining a generally cylindrical scaxmer bore 14 inside of which an 

15 associated imaging subject 16 is disposed. Main magnetic field coils 20 are disposed inside 
the housing 12. The main magnetic field coils 20 are shown simplified to a generally 
solenoidal configuration to produce a main Bo magnetic field directed along a z-direction 
lying parallel to a central axis 22 of the scanner bore 14. The main magnetic field coils 20 
are typically superconducting coils disposed inside in cryoshrouding 24, although resistive 

20 main magnets can also be used. Particularly at higher field strengths and therefore higher 
firequencies, superconducting magnets are preferred. 

The housing 12 also houses or supports magnetic field gradient coils 30 for 
selectively producing magnetic field gradients along the z-direction, along in-plane 
directions transverse to the z-direction (such as along Cartesian x- and y-directions), or 

25 along other selected directions. The housing 12 also houses or supports radio frequency 
head or body coils 32 for selectively exciting and/or detecting magnetic resonances. 
Although birdcage coils are common at 128 MHz and below, other coils besides a birdcage 
coil can be used as a volume transmit coil, such as a transverse electromagnetic (TEM) 
coil, a phased coil array, or other type of radio firequency coil. The housing 12 typically 

3 0 includes a cosmetic inner liner 36 defining the scanner bore 14. Instead of or in addition to 
the radio fi-equency coil 32, a local radio frequency transmit and/or transmit/receive coil 
such as an illustrated TEM head coil 40 can be employed 
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The main magnetic field coils 20 produce a main magnetic field Bo in a 
z-direction which is preferably at least 3.0T, and more preferably greater than 3.0T, such as 
7.0T or higher. A magnetic resonance imaging controller 44 operates magnet controllers 46 
to selectively energize the magnetic field gradient coils 30 and operates a radio frequency 
5 transmitter 50 coupled to one or more of the radio firequency coils 32, 40 to selectively 
energize the radio firequency coil or coils 32, 40. By selectively operating the magnetic 
field gradient coils 30 and the one or more radio frequency coils 32, 40, magnetic 
resonance is generated and spatially encoded in at least a portion of a selected region of 
interest of the imaging subject 16. The magnetic resonance imaging controller 44 operates 

10 a radio frequency receiver 52 coupled to one or more of the radio frequency coils 32, 40 to 
receive magnetic resonance k-space data samples that are stored in a k-space memory 56, 

A reconstruction processor 58 applies a suitable reconstruction algorithm 
such as a Fourier transform reconstruction algorithm to reconstruct the k-space samples 
into a reconstructed image including at least a portion of the region of interest of the 

1 5 imaging subject The reconstructed image is stored in an image memory 60, displayed on a 
user interface 62, stored in non-volatile memory, transmitted over a local intranet or the 
Internet, or otherwise viewed, stored, manipulated, or so forth. The user interface 62 can 
also enable a radiologist, technician, or other operator of the magnetic resonance imaging 
scanner 10 to communicate with the magnetic resonance imaging controller 44 to select, 

2 0 modify, and execute magnetic resonance imaging sequences. 

In one example imaging sequence, the TEM head coil 40 is energized by a 
radio frequency excitation to excite magnetic resonances of in the region of interest of 
the imaging subject 16. For an example main magnetic field Bo=7.0T, the corresponding 
magnetic resonance frequency is 5^es=7Bo=298 MHz, where ^^42.58 MHz/T is a proton 

25 gyrometric ratio. At this frequency, resonance signals have a wavelength in air of about 
one meter, but in the human head with an average permittivity of 8r=40 the wavelength is 
A;^il6cm, which is smaller than many typical imaging regions of interest. During the radio 
frequency excitation, the magnetic field gradient coils 30 produce a slice-selective 
magnetic field gradient along the z-direction to select an axial slice or slab 66 that meets 

30 the resonance condition for excitation. Rather than selecting an axial slice, a coronal, 
sagittal, or otherwise-oriented slice can be selected using a suitably directed magnetic field 
gradient applied during the radio frequency excitation. The resonance signals from the 
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selected slice are read out using a sequence of phase encode gradients applied along an 
in-slice phase encode direction, and readout gradients applied along an in-slice readaut 
direction transverse to the phase-encode direction. In another suitable readout, in-slice 
magnetic field gradients are applied along selected angular directions transverse to the 
5 z-direction and spanning 180° or 360° of angular views. The slice or slab imaging sequence 
is repeated for a plurality of adjacent or spaced-apart slices to produce imaging data that is 
suitably reconstructed by the reconstruction processor 58 into volumetric or multi-slice 
reconstructed image data. The described imaging is only an example. Those skilled in ttie 
art can readily modify the described multi-slice sequence or construct other multi-slice 
1 0 imaging sequences that include features such as spin refocusing, an echo planar readoixt, 
contrast enhancement or selection mechanisms such as inversion recovery pre-pul&es or thie 
like, and so forth. 

With reference to FIGURES 2 A and 2B, the illustrated TEM head coil 40 
includes a plurality of rungs 70. The rungs are shown on-end in phantom in FIGURE 

15 Each rung 70 includes a copper tube 72 having a copper rod 74 inserted therein. In th.e 
illustrated TEM head coil 40, the position of each inserted copper rod 74 is adjustable 
within the corresponding copper tube 72 to tune the resonance fi*equency of the TEM head- 
coil 40. A generally cylindrical copper TEM radio frequency screen 76 surrounds the rungs 
70, and a copper end cap 78 is disposed at one end of the TEM head coil 40. In FIGURE 2, 

2 0 Cartesian x-y-z coordinates are indicated corresponding to the spatial orientation of the 
TEM head coil 40 in the cylindrical scanner bore 14 shown in FIGURE L It will be 
appreciated that the TEM head coil 40 can be used in other orientations; however, only the 
Bi field transverse to the Bo field generates signal. 

FIGURES 3A, 3B, and 3C illustrate a finite difference time domain model 

2 5 for use in electromagnetic modeling of the TEM head coil 40 and the region of interest of* 
the imaging subject 16 which includes the head and upper shoulders. The TEM head coil 
40 is modeled by a TEM head coil model 40' illustrated in FIGURES 3A, 3B, and 3C, 
while the region of interest of the imaging subject 12 is modeled as a head model 16' 
illustrated in FIGURES 3 A, 3B, and 3C. The human head model is based on a human body^ 

30 model that identifies distribution of tissues in the body. In FIGURES 3A, 3B, and 3C, 
various tissues of the human head are diagrammatically indicated by different shadings 
corresponding to differences in conductivity and permittivity characteristics. 
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Hmnan body models and selected animal body models suitable for use in 
finite difference time domain electromagnetic modeling are available from a number of 
sources, such as from: 

5 The United States Air Force Research Laboratory 

Brooks Air Force Base, TX 
http://www.brooks.af mil/AFRL/HED/hedr/ 

and from: 

10 

The Visible Human Project 
National Library of Medicine 
Bethesda, MD 

http://www.nlm.nih.gov/research/visible/visiblejiuman.html. 

15 

Suitable finite difference time domain electromagnetic modeling software can be obtained 
from some of these sources as well, or can be obtained as stand-alone modeling software 
packages into which the human body model 16' and the radio frequency coil model 40' are 
imported. Typically, finite difference time domain electromagnetic modeling operates on 

20 digital models such as pixel-based or voxel-based models of the imaging subject and the 
radio frequency coil. As an example, the human body models presently available from the 
United States Air Force Research Laboratory are voxel-based models having a (1 mm)^ 
resolution. Those skilled in the art can measure and construct other models of human 
beings or other imaging subjects suitable for finite difference time domain electromagnetic 

25 modeling, and/or can use other numerical electromagnetic modeling algorithms besides 
finite difference time domain algorithms for the electromagnetic modeling. Suitable 
conductivity and dielectric data for various tissues of the human body at 298 MHz 
frequency (7.0T main Bo magnetic field) are provided in Table L The source for values in 
Table I is the United States Air Force website 

3 0 http://www.brooks.af mil/AFRL/HED/hedr/reports/dielectric. 



Tissue Type 


o (S/m) 


fir 


p (kgW) 


Blood 


1.690 


58.442 


1058 


Blood Vessel 


0.537 


48.357 


1040 


Body Fluid 


1.518 


69.017 


1010 


Bone (Cancellous) 


0.205 


22.256 


1920 


Bone (Cortical) 


0.096 


14.103 


1990 


Bone Marrow 


0.031 


5.354 


1040 
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Cartilage 


0.552 


46.808 


1097 


Cerebral Spinal Fluid 


2.300 


69.215 


1007.2 


Eye (Cornea) 


1.025 


56.465 


1076 


Eye (Lens) 


0.495 


46.910 


1053 


Eye (Sclera) 


0.955 


59.814 


1026 


Eye (Aqueous Humor) 


1.518 


69.017 


1008.9 


Fat 


0.042 


5.221 


916 


Glands 


0.851 


62.472 


1050 


Ligaments 


0.700 


53.045 


1220 


Lung Inflated 


0.336 


22.645 


260 


Lung Deflated 


0.675 


53.645 


1050 


Lymph 


0.851 


62.472 


1040 


Mucous Membrane 


0.630 


51.957 


1040 


Muscle 


0.916 


66.547 


1046.85 


Nerve 


0.418 


36.951 


1038 


Brain Grey Matter 


0.824 


58.397 


1038 


Brain White Matter 


0.487 


42.871 


1038 


Cerebellum 


0.862 


56.052 


1038 


Skin(I>ry) 


0.534 


43.797 


1125 


Tooth 


0.082 


13.446 


2160 


Model Average 


0.551 


40.058 


1110.02 



Table I - Tissue conductivity and permittivity values at 298 MHz, and density for the 

human head model. 



5 With reference returning to FIGURE 1, the human head model 16* 

including tissue conductivity and dielectric values such as are given in Table I, thes radio 
frequency coil model 40', and an input radio frequency pulse model 80 are input into an 
electromagnetic simulator 82 which implements the finite difference time domain 
electromagnetic modeling or other suitable electromagnetic modeling algoritbm to 

10 compute the spatially dependent Bi field 84 as a fimction of spatial position in two- or 
three-dimensions. For example, the Bi(x,y,z) field is suitably calculated. The spatially 
dependent Bi field computed by the electromagnetic simulator 82 is input to a radio 
frequency transmitter calibration processor 88 which computes, for example, the a^^^erage 
|Bi| value of the Bi field across a slice to generate a per-slice radio frequency transmitter 

15 calibration 92. 

With reference to FIGURES 4 and 5, radio frequency transmitter calibration 
data for the loaded TEM head coil 40 calculated using the TEM head coil model 40' and 
the human head model 16' are shown. FIGURE 4 shows the average |Bi| field in microtesla 
averaged over a slice as a fiinction of axial slice position z. FIGURE 5 shows the 
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normalized |Bi| field standard deviation in naicrotesla from a slice as a function of axial 
slice position z. Also shown in FIGURES 4 and 5 are the corresponding data for 
electrontiagnetic modeling of the loaded TEM head coil 40 using the TEM head coil model 
40* and a homogeneous head model that uses an "average" tissue conductivity and 
5 permittivity for the entire head. In FIGURES 1, 4, and 5, axial slices along the z-dhrection, 
such as the illustrated slice 66, are considered; however, substantially similar processing 
can be employed for per-slice calibration of coronal, sagittal, or otherwise-oriented slices. 

Preferably, the radio frequency transmitter calibration processor 88 also 
conoputes the average, position-dependent, and peak specific absorption rate (S AR) values 
10 from the simulated spatially dependent Bi field 84. Table 11 compares average SAR (in 
units of W/kg) and peak SAR (that is, the maximum SAR averaged over 1-gram portions 
of tissue, or averaged over 10-gram portions of tissue) using the non-homogeneous, 
realistic human head model 16' and using the homogeneous head model that uses an 
"average" tissue conductivity and permittivity for the entire head. 

• 15 



100% 
RF duty-cycle 


SAR average 
over whole 
head model 
(W/kg) 


SAR average 
over the head 
portion inside 
coil (W/kg) 


Maximum 
SAR per 

1-gram tissue 
(W/kg) 


Masimum 
SA.R per 
10-gram tissue 
(W/kg) 


Detailed head 

model 


89.4 


129.5 


831.3 


415.0 


Homogeneous 
head model 


100.1 


148.1 


404.1 


293.5 



Table II - Specific absorption rate (SAR) comparison for non-homogeneous human 
head model and homogeneous head model with same geometry. SAR is scaled to have 
average |Bi|=10 |xT in the central transverse slice of the head model. 

20 

Table II shows one advantage of the non-homogeneous, realistic human 
head model 16' over the homogenous model: the simulated peak SAR values are more 
accurately calculated using the non-homogeneous model 16', since they account for local 
changes in SAR due to different tissue types, tissue masses, interfaces between tissue 
2 5 types, and so forth. In particular, the realistic human head model 16' identifies a factor of 
two larger peak SAR value for 1-gram tissue voxels compared with analogous modeling 
using the homogeneous head model. 
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In one approach for calibrating radio frequency transmissions, the radio 
frequency transmitter caUbration processor 88 computes average and peak SAR values and 
adjusts the radio frequency calibration 92 or advises the user to adjust imaging sequence 
parameters such as echo time, flip angle, repeat time, or number of slices, to ensxire that the 
5 computed average and peak SAR values satisfy relevant regulatory SAR limits 96 relating 
to maximum average and/or peak SAR values. For example, in the United States, the Food 
and Drug Administration (FDA) issues SAR safety regulations, while in much of Europe 
the International Electrotechnical Commission (TEC) regulations are applicable. 

In a suitable approach which substantially maximizes the average field |Bi| 

1 0 while satisfying the SAR regulatory limits 96, a target value for the per-slice average field 
IBil, or an associated parameter such as the associated flip angle, is selected. The radio 
frequency transmitter calibration processor 88 determines a suitable radio frequency 
excitation intensity for each slice that produces the desired average field |Bi| in that slice. 
The calibration processor 88 also computes average SAR, peak SAR, or another relevant 

15 SAR value over the whole excited volume. If one or more of the computed SAR values 
exceeds one of the relevant SAR regulatory limits 96, the target value for the per-slice 
average field |Bi| is reduced, or another imaging parameter affecting SAR is adjusted, and 
the computation is repeated until the SAR regulatory limits 96 are satisfied. This approach 
provides a correct average field |Bi| across the slices while satisfying the safety limits 

2 0 imposed by the SAR regulatory limits 96. 

During imaging, the magnetic resonance imaging controller 44 inputs a 
slice position 100 of the slice to be imaged, such as a z-axis position in the case of axial 
slices, to the radio frequency calibration 92, which then returns the radio frequency 
excitation intensity suitable for that slice. The magnetic resonance imaging controller 44 
25 controls the radio frequency transmitter 50 to apply that excitation to the head coil 40 to 
excite or manipulate magnetic resonance in the imaging subject 12. The radio frequency 
calibration 92 is suitably arranged as a table of radio frequency excitations versus slice 
position, as an empirical fitted mathematical function, or the like derived from the 
simulated spatially dependent Bi field 84. 

3 0 Rather than using the electromagnetic simulator 82 to compute the spatially 

dependent Bi field, a Bi mapping can be acquired experimentally with the region of 
interest of the imaging subject 12 arranged in the head coil 40. A magnetic field map of 
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each slice to be imaged is suitably generated by acquiring and reconstructing a (0, 2G) pair 
of spin-echo images of the slice. A mapping processor 110 performs Bi field mapping 
based on the (G, 29) pair of spin-echo images to generate a per-slice average |Bi| v. z 
calibration 112. This may be performed for a number of subjects to derive a statistica.1 
5 average set of spatial Bi maps for the human head. For imaging using sagittal, coronal, or 
otherwise-oriented slices, analogous suitable per-slice average |Bi| calibrations aire 
similarly generated. During imaging, the magnetic resonance imaging controller 44 
accesses the per-slice average |Bil calibration 112 to determine a suitable radio frequency 
excitation intensity, and the magnetic resonance imaging controller 44 controls the radio 

1 0 fi-equency transmitter 50 to apply that excitation to the head coil 40. 

Rather than performing the Bi mapping on a human head. The Bi map can 
be measured using a spatially non-uniform compartmentalized phantom model of at least 
the region of interest of the imaging subject disposed in the radio frequency coil. For 
example, a spatially non-uniform compartmentalized phantom model of a human head can 

15 be used to generate the per-slice average |Bi[ v. z cahbration 112 for imaging a human 
head. The calibration 112 is suitably stored in tabular form and accessed by the magnetic 
resonance imaging controller 44 during imaging of an actual himian head. Moreover, 
instead of using a Bi mapping, a rapid magnetic resonance tip angle calibration can be 
acquired and used for the radio frequency transmitter calibration. A combination of these 

2 0 methods may also be used. The rapid tip angle calibration may be used to determine an 
initial calibration relationship between RF transmitter input level and actual average Bi 
value achieved at a given slice position. This data set provides a reference level from 
which radio frequency transmitter input levels of the remaining slices are calculated based 
on the average Bi values modeled or calculated for those slices. 

2 5 The per-slice radio frequency intensity calibration has been described witli 

reference to the illustrated TEM head coil 40 and for imaging using axial slices. However, 
the per-slice radio frequency calibration is readily appUed to substantially any type of radio 
frequency coil, such as a birdcage coil, a phased coil array, a SENSE coil array, or the like, 
and is readily applied to head coils, whole body coils, local coils, and the like. Still further, 

3 0 the described per-slice radio frequency calibration is readily applied to any imaging slice 

orientation, such as a sagittal orientation, a coronal orientation, or a non-standard slice 
orientation. 
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In the illustrated embodiment, the Bi field was characterized on a par-slice 
basis by an average |Bi| value for the slice. That is, the |Bi| field averaged over an area of 
interest of a slice served as the figure of merit for determining the radio frequency 
transmitter calibration 92. However, other figures of merit can be employed, such as a 
5 median |Bi| value, or a root-mean-square |Bi| value given by: 



|l|B.(x,y)r 



Figure of merit = "^j-^^^ ^^^^ (1) 



where in Equation (1), N is the number of slice pixels in the area of iaterest of the slice and 
10 the sunmiation is performed over all N slice pixels. The slice pixels are designated by 

coordinates (x,y), although non-Cartesian coordinates can also be used. 

Other statistical aggregations besides mean, median, average, or 

root-mean-square can be used. The area of interest that is aggregated may include the 

entire field of view of the slice, or the area of interest may include only those pixels 
15 exhibiting a magnetic resonance signal greater than a threshold value indicative of 

corresponding to human tissue, or the area of interest may be an area of pre-selected 

dimensions. 

Although an aggregation figure-of-merit is typically preferred to 
approximately accoimt for in-slice Bi field variations, an unaggregated figure of merit, 
2 0 such as the Bi field along the central axis of the coil with an accounting of the effects of 
the human head, can be used. An unaggregated figure-of-merit may have computational 
efficiency advantages. 

When the volume defined by the areas of interest of the slices over which 
the figure of merit is calculated is less than the total excited volume of the imaging subject, 
2 5 the average and peak SAR values are preferably computed with respect to the total excited 
volume of the imaging subject to ensure that the SAR does not exceed the regulatory SAR 
limits 96 an3nvhere in the imaging subject. 

Moreover, the figure of merit can be the spatially varying Bi field, for 
example, the total range of Bi field amplitude over a slice, and the radio frequency pulse ' 
30 shape 80 can be an adiabatic pulse that is adjusted based on the spatially varying Bi field 
figure of merit such that the adiabatic pulse substantially flattens the in-slice spatial 
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variations in the resulting flip angles. The amplitude and other characteristics of the 
adiabatic pulse are also adjusted such that the substantially flattened flip angles are 
substantially similar for each slice. In this case, the SAR calculation should be performed 
for each adiabatic pulse to ensure that no slice excitation adiabatic pulse will produce an 
5 SAR that exceeds the regulatory limits 96. 

In one embodiment suitable for multi-stage imaging methods such as 
magnetic resonance angiographic peripheral run off studies or whole body 
scanning/screening studies, the imaging subject 16 is moved axially, that is, transverse to 
the acquired slice or slices, continuously or in discrete steps through the bore 14 along the 

1 0 bore axis 22. Imaging is repeatedly performed at a selected bore slice position along the 
bore 14 or over a volume slab of bore slices to perform multi-slice magnetic resonance 
imaging. In this embodiment, the position of the imaged bore slice or slab of bore slices is 
referenced to the bore 14 and does not change; however, as the imaging subject 16 is 
moved axially the portion of the subject 16 coinciding with the selected bore slice or slab 

15 of bore slices does change. 

In such multi-stage imaging methods, the per-slice Bi field calibration can 
be repeated for the selected stationary bore slice each time the imaging subject is moved, to 
account for changes in the Bi field characteristics of the stationary bore slice caused by the 
different portion of the imaging subject coinciding therewith. Each per-slice Bi field 

2 0 calibration can be adjusted based on the previous bore slice calibration, since the change in 
Bi field characteristics of the selected bore slice due to gradual or small-step motion of the 
imaging subject is typically small. In the case of continuous axial subject motion, each 
per-slice Bi calibration computation can be timed so that the calibration is performed for a 
slice of the imaging subject 16 that will coincide with the selected bore slice at the time of 

2 5 the next slice data acquisition. 

The invention has been described with reference to the preferred 
embodiments. Obviously, modifications and alterations will occur to others upon reading 
and understanding the preceding detailed description. It is intended that the invention be 
construed as including all such modifications and alterations insofar as they come within 

3 0 the scope of the appended claims or the equivalents thereof. 
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CLAIMS 

Having described the preferred embodiments, the invention is now claimed to be: 

1. A magnetic resonance imaging method for performing multi-slice magnetic 
resonance imaging of a region of interest of an associated imaging subject (16) using a 
radio frequency coil (40) arranged to generate a Bi magnetic field in the region of interest, 
the method comprising: 

determining a per-slice Bi field value for each slice that is representative of the Bi 
field over a selected area of the slice; 

determining an adjusted per-slice radio frequency excitation for each slice that 
adjusts the Bi field value for the slice to a selected value; 

acquiring magnetic resonance imaging data for each slice using the adjusted 
per-slice radio frequency excitation for that slice; and 

reconstructing the acquired magnetic resonance imaging data into a reconstructed 
image representation. 

2. The method as set forth in claim 1, wherein the determining of a per-slice Bi 
field value for each slice comprises: 

determining the Bi field across at least the selected area of the slice; and 
computing an aggregate value of the detennined Bi field across the selected airea of 
the slice. 

3. The method as set forth in claim 2, wherein the computing of an aggregate value 
comprises: 

computing an average, median, or root-mean-square value of the determined Bi 
field across the selected area of the slice. 

4. The method as set forth in claim 2, wherein the determining of the Bi field 
across at least a selected area of the slice comprises: 

computing the Bi field numerically using a model (40') of the radio frequency coil 
(40) and a non-homogeneous model (16') of the imaging subject (16), the 
non-homogeneous model (16') of the imaging subject (16) eniploying different 
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conductivity and permittivity values for different materials of which the imaging subject 
(16) is formed. 

5. The method as set forth in claim 4, wherein the non-homogeneous model (16*) 
of the imaging subject (16) is a non-homogeneous model (16') of at least a portion of a 
human imaging subject (16), the model (16*) employing different conductivity and 
pennittivity values for different types of tissue. 

6. The method as set forth ia claim 2, wherein the determining of the Bi field 
across at least the selected area of the slice comprises: 

measuring a Bi map of at least the selected area of the slice with one of (1) the 
region of interest of the imaging subject (16) disposed in the radio frequency coil (40), (2) 
the region of interest of a representative distribution of imaging subjects acquired u priori 
in the radio frequency coil (40), and (3) a spatially non-uniform compartmentalized 
phantom model of at least the region of interest of the imaging subject disposed in the 
radio frequency coil (40). 

7. The method as set forth in claim 1, wherein the determining of a per-slice Bi 
field value for each slice comprises: 

determining a value of a figure of merit for each slice that is representative of the 
Bi jSeld over a selected area of the slice. 

8. The method as set forth in claim 7, wherein the determining of an adjusted 
per-slice radio frequency excitation for each slice comprises: 

determining an adjusted radio frequency excitation for each sUce that adjusts the Bi 
field figure of merit to a selected value of the Bi field figure of merit that is substantially 
the same for a plurality of slices for which imaging data is acquired. 

9. The method as set forth in claim 1, wherein the determining of an adjusted 
per-slice radio frequency excitation for each slice includes determining an adjusted radio 
frequency excitation for each slice that adjusts the per-slice Bi field value to a selected 
value that is substantially the same for a plurality of slices for which imaging data is 
acquired, the method further comprising: 
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detemiining a specific absorption rate based on the adjusted per-slice radio 
frequency excitations; and 

conditional upon the determined specific absorption rate exceeding a regulatory 
safety limit, repeating the determining of an adjusted per-slice radio frequency excitation 
for each slice using one of a lower selected value of the per-slice Bi field value and 
adjustment of at least one other imaging sequence parameter to reduce the specific 
absorption rate. 

10. The method as set forth in claim 1, wherein the adjusted per-slice radio 
frequency excitations are adiabatic radio frequency excitations, and the determining of the 
adjusted adiabatic radio frequency excitations comprise: 

for each slice, computing an adjusted adiabatic radio frequency excitation that 
substantially corrects for a variation of the Bi field across the selected area of the slice to 
provide more uniform flip angles. 

11. The method as set forth in claim 1, ftirther comprising: 

moving the associated imaging subject (16) in a direction transverse to the slices, 
the determining of a per-slice Bi field value, the determining of an adjusted per-slice radio 
frequency excitation, and the acquiring of magnetic resonance imaging data being repeated 
for a stationary slice position with the imaging subject (16) moved relative to the stationary 
slice position between each repetition. 

12. The method as set forth in claim 11, wherein the moving of the associated 
imaging subject (16) is one of: 

in discrete steps, with the subject (16) motionless during each 
repetition of the acquiring, and 

continuous, with the subject (16) moving during each repetition of 
the acquiring. 

13. A magnetic resonance imaging apparatus comprising: 

a main magnetic field coil (20) generating a main magnetic field; 

magnetic field gradient coils (30) selectively generating magnetic field gradients; 
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a radio frequency coil (32, 40) arranged to generate a Bi magnetic field in a region 
of interest of an associated imaging subject (16); 

a radio frequency transmitter (50) selectively energizing the radio frequency coil 
(32,40); 

a radio frequency receiver (52) selectively sampling the radio frequency coil (32, 
40); and 

a processor progranamed to perform the method of claim 1. 

14. A magnetic resonance imaging apparatus for performing multi-slice magnetic 
resonance imaging of a region of interest of an associated imaging subject (16), the 
apparatus comprising: 

a radio frequency coil (40) arranged to generate a Bi magnetic field in the region of 
interest; 

a means (82, 88, 110) for determining a per-slice Bi field value for each slice that is 
representative of the Bi field over a selected area of the slice; 

a means (92, 112) for determining an adjusted per-slice radio frequency excitation 
for each slice that adjusts the Bi field value for the slice to a selected value; 

a means (10, 44, 46, 50, 52) for acquiring magnetic resonance imaging data for 
each slice using the adjusted per-slice radio frequency excitation for that slice; and 

a means (58) for reconstructing the acquired magnetic resonance imaging data into 
a reconstructed image representation. 

15. The apparatus as set forth in claim 14, wherein the means (82, 88, 110) for 
determining a per-sUce Bi field value for each slice comprises: 

an electromagnetic simulator (82) receiving a digital model (16') of the region of 
interest and the digital model (40') of the radio frequency coil (40) and estimating the Bi 
field generated across the region of interest, the digital model (16') of the region of interest 
mimicking non-imiform dielectric and conductivity properties of the region of interest. 

16. The apparatus as set forth in claim 15, wherein the electromagnetic simulator 
(82) employs a finite difference time domain algorithm. 
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17. The apparatus as set forth in claim 14, wherein the means (82, 88, 110) for 
determining a per-slice Bi field value for each slice comprises: 

a means (82, 88, 110) for determining a value of a figure of merit for each slice that 
is representative of the Bi field over a selected area of the slice. 

18. The apparatus as set forth in claim 17, wherein the means (82, 88, 110) for 
determining a value of a per-slice Bi field figure of merit for each slice that is 
representative of the Bi field over a selected area of the slice comprises: 

a statistical aggregation means (88) for calculating an aggregate value 
representative of the Bi field over the selected area of the slice. 

19. The apparatus as set forth in claim 14, wherein the means (92, 112) for 
determining an adjusted per-slice radio fi*equency excitation for each slice comprises: 

a table (92, 112) of adjusted per-slice radio frequency excitation versus slice 
position derived from the per-slice Bi field values. 

20. The apparatus as set forth in claim 14, further comprising: 

a means (88) for determining a specific absorption rate based on the adjusted 
per-slice radio firequency excitations; and 

conditional upon the determined specific absorption ratio exceeding a regulatory 
safety limit, repeating the determining of an adjusted per-slice radio frequency excitation 
for each slice using one of lower selected values of the per-slice Bi field value and 
adjustment of at least one other imaging sequence parameter to reduce the specific 
absorption rate. 
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